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Abstract 

We add massive fundamental matter to the ABJM model by adding 
D6-branes wrapped asymptotically over RP^. We find two types of 
solutions at finite temperature, one that enters the black hole and one 
that ends before the black hole. We analyze the behavior of the free 
energy as a function of temperature, and find that the system exhibits 
a phase transition between the two types of solutions, similar to what 
happens in the D3-D7 system. We also analyze the meson spectrum 
in the model and find several massive scalar modes, again, quite like 
the D3-D7 system. We end with a calculation of the conductivities in 
the two phases. 



1 Introduction 

The ABJM |l] model is a conjectured holographic duality between a three 
dimensional A/" = 6 supersymmetric Chern-Simons theory with gauge group 
U{N)k X U{N)^k and M-theory in an 74^5*4 x S"^ /Z^ background or type 
IIA string theory on AdS/i, x CP^ . Holographic dualities have been useful in 
understanding strongly interacting systems. Most of the applications have 
been to four dimensional systems, such as N=4 SYM. The ABJM model 
provides an opportunity to study strongly interacting three dimensional sys- 
tems, which might be applicable to phenomena in condensed matter physics. 
It is therefore important to extend the model so that it can represent more 
realistic systems containing matter degrees of freedom. 

The basic ABJM model contains non-dynamical gauge fields (since it is 
a CS theory) and several scalar and fermionic fields in the bi-fundamental 
representation. There is no field which is solely in the fundamental of only 
one of the groups. 

It is possible to extend the model to include fundamental matter by 
adding D6-branes that wrap an RP^ subspace of CP^ in the Type IIA de- 
scription [2]-[5]. In this paper we will study D6-brane embeddings corre- 
sponding to massive matter and explore the behavior of the system at finite 
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temperature. The D6-brane is treated as a probe in the AdS-Schwarzschild 
Type IIA background. Massive embeddings have been considered previously 
at zero temperature by Jensen [6] as a model for a BKT type transition in 
a background magnetic field. I will show that at finite temperature there 
are two kinds of embeddings, corresponding to two possible phases of the 
dual field theory. In a black-hole (BH) embedding the D6-brane intersects 
the horizon, and this corresponds to an ungapped phase in the field theory. 
In a Minkowski (MN) embedding the D6-brane ends outside the black-hole 
horizon, and this is dual to a gapped phase. The MN embedding is preferred 
at low temperature. As the temperature increases there is a first-order phase 
transition to a BH embedding. 

The behavior of the ABJM-D6 system is qualitatively similar to the D3- 
D7 system studied in [7]-[9], and can be viewed as the three-dimensional 
version of that system. 

The rest of the paper is organized as follows. In section 2 we will briefiy 
review the ABJM model and explain how to extend it to include matter 
in the fundamental representation. In section 3 we will derive and solve 
the equations of motion, and then in section 4 we will exhibit the phase 
transition. In section 5 we study the meson spectrum. In section 6 we 
analyze the electrical conductivities. We end in section 7 with some short 
conclusions. 

2 The geometry of the model 
2.1 The Type IIA Background 

The ABJM model describes M2-branes on a C^/Zk orbifold singularity. 
The low energy field theory is three dimensional A/" = 6 supersymmetric 
Chern-Simons gauge theory with a gauge group U{N)k x U{N)^k- The dy- 
namical fields are four massless complex scalars in the bi-fundamantal rep- 
resentation and their fermionic partners. At large A^ this theory has a dual 
supergravity description corresponding to the near-horizon background of the 
M2-branes. For k <^ N^^^ this is given by eleven dimensional supergravity 
in AdSi X S'^ jZ^,. For A^^/^ <^ k <ti N the more appropriate description is in 
terms of Type IIA supergravity in AdS4^ x CP^. This will be the regime in 
which we work. 
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The AdS^ x CP3 metric is given by: 

= l^^Ads, + dslp^ \ (1) 
ds\ds, = ^ [-dt" + dx^ + dy^) + ^dr^ , (2) 
where R is the radius of the S"^ in lid Planck units and is given by [l]: 

R^ = 2hVNk. (3) 
In a standard parameterization the CP^ metric is given by: 

d^cPs = dc, + cos 4 sm c,{dip H ^ — a0i ^ — ct02j 



where the ranges of the angles are as follows: 

< ^ < 7r/2 , < ^ < 27r , < < vr , < (f)i < 27t . 
The dilaton and RR fields are given by [l]: 

= % (5) 
F2 = kJ (6) 



3 3 /r\5 

F4 = -R^e^ = -R^ - dtAdxiA dx2 A rfr (7) 
8 8 V 2 / 



3„, 3„, /r\5 
-2. 

where J is the Kahler form of CP^, and is given explicitly by 



J = — cos ^ sin A (2(i'0 + cos — cos ^2*^02) 

1 o 1 



cos^ ^ sin 6'i(i6'i A dcpi — - sin^ ^ sin 92d92 A 



2 " ^ ^ 2 

There is a generalization of the ABJM model (known as the ABJ model), 
which allows for a diffrence, /, between the ranks of the two groups flO]. In 
this case the Type IIA dual has a non-zero B-field given by: 

B2 = [j. (9) 



3 



As we are interested in the behavior of the system at finite temperature, 
we will need the metric when temperature is taken into account. At finite 
temperature the relevant background contains a black hole, and the metric 
becomes: 



where rj^^ is the horizon radius given by: r^j 



I&ttT 



2.2 Adding flavor 

Now we shall describe the embedding. First, we should mention that once we 
reduce from M-theory to type IIA the M2 branes become D2 branes living 
on the t ,Xi , X2 coordinates of AdS^. We want to add fundamental matter 
by adding a single probe D6-brane, so that the strings between it and the D2 
branes will give us fundamental fields. The distance between these branes 
will determine the mass of these fields. 

Our starting point is the massless embedding (2]-[5], in which a D6-brane 
is embedded on all of AdS^^ and on an RP^ subspace of CP^, which is the 
Lagrangian submanifold. The RP^ d CP^ that the D6-brane wraps is given 
by: 

6^ = 6^ = 6, 01 = -02 = 0, e = ^- 

This embedding preserves TV = 3 supersymmetry. Since the D6-brane shares 
all the coordinates of the D2-branes, the matter fields live in all 3 dimensions 
of the dual field theory. This gives the field theory with the quiver diagram 
shown in figure 1. 

We will consider a more general embedding given by ^ = ^{r), with 
^{r — 7- 00) = |. This corresponds to a mass deformation for the matter 
fields. In order to compute the mass we need to undo the large limit and 
measure the D2-D6 separation. This is not easy to do directly in the Type 
IIA description, since the seven-dimensional geometry transverse to the D2- 
branes at finite A^ is not well-understood. However, we can start with the fiat 
space Type IIB brane configuration of [l] and work our way to the Type IIA 
picture. Recall that the Type IIB configuration consists of an NS5-brane, a 
(l,fc) 5-brane and D3-branes arranged as follows: 
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Figure 1: Quiver diagram for the field theory with flavor. 



0123 4 567 8 9 
NS5 

(1, fc)5 • • • cos6 cos 6 cos 6' sin^ sin^ sin^ 
D3 • • • • 

In this picture fundamental matter is incorporated by adding a D5-brane 
along (012789) [2]-[5]. Since the background spacetime is flat in this descrip- 
tion, the mass of the fundamental is given simply by the D3-D5 separation 
along (345). The Type IIA supergravity dual is obtained by first T-dualizing 
along x^, then lifting to M theory, taking the near- horizon limit, and finally 
reducing along the Hopf fiber to Type IIA string theory. The transformation 
of coordinates is given explicitly by [2]: 

x'^ = ^ (11) 
x[ = Xi = r"^ cos^ ^(cos6'i, sin 6*1 COS01, sin 6^1 sin 0i) (12) 

X2 = Xi + kx2 = r'^ sin^ ^(cos 62, sin 6*2 cos 02, sin 6*2 sin02) (13) 

where xi = {x'^,x^,x^) and X2 = {x^,x'^,x^). Furthermore, since the above 
relation involves T-duality transverse to NS5-branes, there is a non-trivial 
metric given by [l]: 



^ 2 doc ' doc ^ doC'2 ' dx ^X4^ 
2 1 3^ -j^ I 21^2! 

In this description the mass is given by the proper length of a (geodesic) path 
of fixed xi, i.e., 
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^ - f ^7= - V2{.M+k^\ - (15) 
J V2F2I 

= "\/2r^(sin^ - cos^) = 2r^sin(^ - ^) 
The mass is given by 

M^'±^. (16) 

This embedding preserves some of the symmetries of CP^. The asymp- 
totic space, i?P^, has an 5'0(4) symmetry which is locally equivalent to 
5*0(3) X 50(3). This embedding preserves one of the 50 (3) 's (which corre- 
sponds to the R-symmetry in the field theory) and breaks the other to U (1) 
(which is a global symmetry in the field theory). 



3 Derivation and solution of the equation of 
motion 

Next we determine the behavior of ^(r) from the equation of motion of the 
brane. The relevant part of the D6-brane action is given by Sdbi + Scsi 
where: 

Sdbi = -A*6 j e~^^/-det{Gab)dJo^ (17) 

5c5 = -/i6 J P[Cr]. (18) 

The 7-form potential can be evaluated by dualizing F2 to get Fg and then 
integrating to get O7. The process is not difficult but quite tedious so we will 
simply write the result: 

dCj = -———J A J Adr Adt Adxi Adx2. (19) 
In a specific gauge choice: 

R^r^ 

Cr^ -^-^^J AJ AdtAdXiAdx2, (20) 
giving the following puUback to the brane worldvolume: 

R^r^ 

P\C^] = -— ——cos 2^ sin 2^ sin A d6l A (21) 
3 • 2^^A;^ 
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The induced metric on the D6-brane is: 



,2 _ ^ 



r 

16 



+d9'^ + sin^ 9d(f)'^ + sin^ 2^{dij + cos 
The D6-brane worldvolume action then reduces to: 

3=^1 drr' sm{20 ( v^r-^ + - ^ cos(20e') 
and the EOM for ^ can be easily derived: 



'2 » dr' 
>)'] -(22) 

(23) 



5r 



r^sin(20f (1 



r.6 , 



l + e^r2(l-^) 



2r^cos(20 Wl + e^r2 1 



H 



sin(2e) . 



(24) 



Note the trivial solution ^ = 7r/4, corresponding to the massless embedding. 
The asymptotic form of the general solution is found to be: 



TT 



(25) 



with «_ = —4 and = — 2j^ The leading term is proportional to the mass, 
where the constant can be calculated using the mass formula (16). The 



normalizable mode, in accordance with the state-operator correspondence, is 
related to the expectation value of the dual operator. Therefore, this term 
corresponds to the flavor condensate c. 



From (22) one can see that the brane can end at some ro > rjy, where 
or |. Indeed, as was mentioned, there are two possible cases: 
either the brane reaches the horizon or it ends before it. The behavior of the 
solution in each case is different. 

The equation of motion can be solved numerically for both types of solu- 
tions. In the black hole embedding, one starts from the horizon with some 
initial value ^{th), which determines the mass. The derivative at the horizon 
^'(rn) is fixed by the EOM: 



^As in other AdS-CFT examples, this is related to the dimension of the dual operator. 
Note however that our radial coordinate is different than the one related to energy scaling. 
The relation being VEnergy — {fHereY , so the dimensions of the field are 1, 2 consistent 
with mass^ = — 2 and with the results of previous works 
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cot 2^0(1 + sin 2^0 



3r 



H 



(26) 



In the Minkowski embedding one starts at some r = ro > rn, which now 
determines the mass, and ^'(ro) — ?• 00 (in practice we use vanishing derivative 



of r(^) ). Once the solution is determined one calculates the mass from (25). 
The results we found are shown in figure 2. 




(a) 




Figure 2: ^ as a function of initial parameter at fixed temperature: (a) 
BH solution as function of the initial angle, (b) MN solution as function of 
2) 



the initial r. {r^ 



Since the equations are invariant under ^ — )■ | — ^, every embedding (MN 
or BH) has a "mirror" embedding, but with the opposite sign of the mass. 
This is apparent in figure la, which shows the BH solutions. Figure lb shows 
MN solutions with ^(ro) = 0. There are corresponding ones for ^{ro) = |. 



4 Phases of the model 

Looking at figure 2, we see that when M » T there is only one solu- 
tion, which is an MN embedding far from the horizon. As the temperature 
increases two BH solutions and another MN solution appear. When the tem- 
perature is further raised the MN solutions disappear and we are left with 
only a single black hole solution. This suggests that there should be a phase 
transition between these two states at some temperature. We next set out 
to confirm this. One way to do this is to examine the c vs. ^ curve which 
we have calculated. The calculation is quite similar to the evaluation of the 
mass range, but now we used the value of ^ and its derivative at large r to 
calculate M and c by solving the following two equations: 
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^(r oo) 



^'(r oo) 



71 

4 



ttM c 
+ -r 



ttM 



4- 



(27) 



(28) 



The results are shown in figure 3. 

First, one notices again the ^ — )■ | — ^ refiection symmetry where on the 
right are embeddings with ^ < | while on the left ^ > f • More importantly, 
there is a region where the graph is not single-valued, this is typical in cases 
of a first order phase transition, and confirms that this also happens in this 
case. Furthermore, one can use the Maxwell construction to calculate the 
critical value of ^. 
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Figure 3: c vs y" constant temperature: blue and green are MN embed- 
dings while the purple are the BH embeddings. The Maxwell construction is 
illustrated so one can infer the critical value of 

Another way to show the phase transition is to evaluate the free energy 
of each solution, for a fixed mass, as a function of temperature. The free 
energy of a solution is proportional to its Euclidean action: 



F oc Se- 



(29) 



This diverges at large r. Therefore, we must regularize it. We used holo- 
graphic renormalization, adding the following counterterm: 



S, 



counterterm 



I d'xv^ (i - 2(e - J?) 

J bound. '-t / 



(30) 
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where 7 is the induced metric on the UV cutoff. This renders the action 
finite. However since there is a CS term one must also add a boundary term 



so that the action will be gauge invariant. Once the gauge choice (20) is 
made there remains the gauge symmetry: 

= "3T2np^°^2esin2a' ^ ' 3 . 211^2 ^^^^ 
Hence, we add the following boundry term: 

Shoundary = I d? X — SIV? {2^) , (32) 

^ Jbound. ^'^ 



SO that under (31) it transforms as: 



^/ d».^sm^(2f)^^/ d».!:^sm^(20. (33) 

^ Jbound. -L^ Z ft. Jbound. 

and the action remains invariant. We can now use this gauge symmetry to 
set C^{r = th) = 0, and so we have: 

= -^^(1 - ^)cos2esin2ee', (34) 

Sboundary = ^'^^(1 " ■ (35) 

^ Jbound. ' 

We then scanned over a range of temperatures, evaluated the solutions 
with the appropriate mass for each temperature and calculated their free 
energy. A typical plot of the free energies is shown in figure 4. 

As expected, one can see that at low temperatures the MN solution is the 
only possible one. As the temperature increases, two BH solutions appear, 
but are unfavored. Then at some temperature one of the BH solutions be- 
comes favored and a phase transition occurs. Also, at some temperature the 
other BH solution pinches off and becomes an MN solution (one can see in 
figure 2 that the masses for ^{th) ~ in (a) and tq ^ th in (b) are nearly 
identical), which eventually merges with the original MN solution. At high 
enough temperatures there is only the BH solution. 



5 Meson spectrum 

We will now analyze the meson spectrum corresponding to the fluctuations 
of the probe brane worldvolume fields. There are three scalar fields related 
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Figure 4: Free energy vs temperature: Blue and orange are MN embeddings 
while red and green are BH embeddings 



to the directions normal to the brane, and there is the gauge field which 
gives one scalar field and one vector meson. Furthermore, there are various 
fermionic fields which are the superpartners of the bosonic fields. I will look 
at only the first three scalar fields, corresponding to the three embedding 
coordinates. We will use the following notation: 

^ = ^0 + 5^, — t; — = — ^ — = (36) 



where is the solution of equation (24), and it is assumed that 6^, a, j3 « 1 



and that they are independent of the CP^ coordinates, but may depend on 
the AdSi coordinate^ Next, we need to evaluate the probe brane action, 
which will now depend on the fields 6^, a, /3. As these fields are very small 
we can expand the action to the first non- vanishing order (which is quadratic 
as we expand around a solution to the probe brane's EOM). We will have a 
system of coupled hnear second order differential equations. The eigenvalues 
of these equations are the masses of the mesons. 

However calculating the full action when all three fields are present is 
quite tedious. This can be eased by noticing that (3 decouples from the 



^Dependence on the coordinates of CP^ will simply give angular momentum internal 
modes for each meson. 
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rest. The reason is that the metric is independent of 0i,02, and therefore /3 
enters the equations only through derivatives. These are either second order 
and diagonal in /3 (like [dP)"^), or first order terms but non-diagonal (like 
dpd(f), which will give a term on AdS4,(j)). The diagonal terms are already 
second order in /3 and therefore cannot couple. The non-diagonal terms 
always multiply one another (because the original AdS^ metric is diagonal) 
generating terms which are second order in /3. This makes f3 easier to analyze 
and so we shall start with this field. 
Therefore we take: 

01-02 <t>\+ (t>2 , . 

e = eo, ^Y^ = 0' ^^^i^ = /3«l. (37) 

Expanding the Lagrangian to second order gives: 



Lo + VLr^ sin^ 2^o 



6 



16 



1 



{31 



Where Lq is the zeroth order lagrangian giving the action (23), and fl is some 
numerical constant. The equation of motion for f3 is then: 



dr 



r ' sin^ 2^0 ( 



1+ 1 



«)2 



64r sin^ 2^0-1/ 1 + 1 



'\2 



1 H 



(39) 



We solve the equation by separation of variables: 



/3 = /3(r)e- 



-iuit A]i--x. 



(40) 



We then get the following equation for /3(r) 



'1+(1- 



)K,)2 



64r sin^ 2^oa/1 



1 - 



Ho)' 



IkP- 



/3. 



(41) 
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We then need to give either a; or |k| and determine the other such that 
the solution is normahzable. Since the system at finite temperature is not 
Lorentz invariant, one must be careful about the definition of the rest mass. 
We follow [9] and define the rest mass as m = in the rest frame, that is 
when |k| =0. Therefore we set |k| = and solve for the values of u for 
which P{r) is normalizable. 



First we evaluate the asymptotic r 00 dependence, by expanding (41) 
to the first non- vanishing order in r~^, finding: 

(3{r ^ 00) ^ A + Br~^. (42) 

There should therefore be one normalizable solution and one non normaliz- 
able solution. 
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0.0030.0040.0050.0060.0070.0080.009 



Figure 5: The mass of the first mode of /3 as a function of the embedding 
mass, for MN embeddings 

The equation is then solved numerically. First, we shall discuss the MN 
embedding. What we did in practice is to integrate the equation starting from 
a small value of /3(r), with vanishing derivative, near r = tq. The solution 
then diverges at large r, going to: (3{r — )■ 00) — j- ±00. We scanned over u 
and identified the normalizable mode by the value of u where the asymptotic 
behavier of (3 changes from /3(r — i- 00) — )■ +00 to /3(r — )• 00) —00. The 
result for the lowest mode of /3 is shown in figure 5. 

In the BH embedding there are no normalizable modes, as can be seen 



by expanding the EOM (41) around r ^ and solving for the behavior of 



(3 in that region. One finds: 
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P ^ A cos In (r — th)^ + B sin In (r — rj^)^ . (43) 

This is ill-defined as r — )> r and therefore the solutions are non-normalizable. 
This makes sense, since in the BH case, there is nothing preventing the 
mesons from falling into the horizon and hence deeming them unstable. 

Next we will analyze ^ and a. Though it isn't obvious a-priori, the two 
fields actually decouple and so, to simplify matters, we will expand them 
separately. We will start with a and take (as /3 also decouples): 

With this choice the determinant can be broken into two pieces: an AdS + 
9 part and a. ip + 4> part. We will start with the ip + 4> part, which gives: 

sin 2f n I 

J ysin2^ + Q;sin2^cos2^ + Q;2cos2^ci^. (45) 

Expanding the square root to quadratic order gives: 

— ^ / sin^ l + a;cot^cos2£+ ^]d9. 46 

4 7o V 2sin2^; ^ ^ 

Note however that the last term gives a divergent contribution. Doing 
the integral first and then expanding does not work either. The integral can 
be expressed as an elliptic integral but cannot be expanded around a = 0. 
There appears to be a problem at = 0, vr. This can be understood from the 
definition of 9, a and from the range of the angles 9i,92. The only way to 
have 6* = is if 6'i = 6*2 = and therefore a — 0, and likewise for 9 = n. So 
9 — 0,7r and a 7^ is inconsistent and this is the reason for the divergence. 
a must depend on 9 such that a = at ^ = 0, tt. This does not change the 
splitting of the determinant or the ip -\- (p psxt we calculated. Now we need 
to calculate the AdS + 9 part, expand it in a and multiply it by the ijj + (j) 
part. This is quite lengthly and uneventful, and in the end we get: 



L[q;] = Qr^sin^ 2^0-4/1+ 1- 



H 



I 1 



# a" sin 9 



.8(1 + (1 - J ) (^^0)^) 



8sin^ 



a 



- {d,,ay - {d,,ay + - sin 9{d0ay + 



a 



2 sine 



cos 29 
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where I have used a dot and a prime for the time and r derivative respectively. 
We can now derive the EOM for a: 



dr 



.4 

^6 



1+ 



= sin^ 2^0 



64rWl 



l\2 



Proceeding by separation of variables we will take a = aA{r,t, Xi, X2)a0{9) , 
and after some rearrangement we get: 



64 



-dr 



"^sin3 2eo(l-J^) a'j, 



a 



1 - 



- ^4" 



4 cos 2^ de{sm.9deae) ,2 



sin^^ 



sin Oae 



k\ (49) 



We get the following two equations for the two parts: 



dr- 



sin3 2^0 



+64r sin=^ 2eo^l+(l- J) {rQ' 



sin' 2^0 W 1 + (1 - J) {rQ^aA 



''h 



(50) 



'cos 29 

ddismOdeaQ) = ( — + k^sinO 1 ag. 

^ smO 



(51) 



As one can notice by inspection, ae — C sin ^ is a solution of the last 
equation, with k"^ = 0, obeying the boundary conditions. It is the lowest 
possible mode of ae, which is what we are after. One can explore higher 
modes but we will not do it here. With this, the first equation becomes: 
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dr 



■2«o 1 



1+1 



«)2 



64r sin^ 2^c 



1 - 



-M- 



, (52) 



which is identical to equation (41). So, we find that a and (3 have identical 



masses at least for the lowest a modes in CP^. This can also be inferred 
from the symmetry of the embeddings. As previously mentioned, part of the 
symmetry of RP^ is an S'0(3) symmetry which becomes a global symmetry 
in the field theory. The three scalars are a fundamental multiplet of this 
symmetry. In the massive embeddings this symmetry breaks to a U{1) = 
SO (2) under which the fundamental multiplet breaks into a singlet and a 
two-vector. The fields /3, a belong to the two-vector and therefore they still 
have the same mass, even to higher order. 

Now we turn our attention to ^. Expanding it as ^ = + doesn't 
work. That is because fiuctuations in the MN embedding need also to take 
into account slight movement in tq, however the boundary conditions forces 
^(ro) = |,0 and therefore cannot allow such movement. So we need to 
convert to a different variable. We define a new radial coordinate p: 



and new "Cartesian" coordinates m, y. 

1 ( ^ 
u = p cos I 24 — — 



y = p sin 2^ 



TT 



(53) 



(54) 



(55) 



We will describe the fluctuations in terms of y{u) \6\. Now we can expand 
the flelds: 



y = yo + 5y, 



01 — 02 



0, 



+ 02 



0. 



(56) 



2 ' 2 

After evaluating the AdS determinant and the Chern-Simons term one 
gets the following Lagrangian: 
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m 




M 



Figure 6: The mass of the first mode of y as a function of the embedding 
mass, for MN embeddings 



L[5y] = Lo + n^ 



(l-pl2)(V)2 
2p8(Vl+{?/')2)3 



(57) 



where 



f{u,y,y') 
and 



yy' I A 



,12 



3p 



12 



6 - ( 1 + 



(5^ 



iu,y,y') = y^T+Wy 



35. 

pi2 p4l +pl2^ 



Qy^juy' - y) 



1 + 
■ 3pi2 



12 - 2mV^ ( 1 + 



3y{uy' - 
p2 



2^, 2^7iu'-2y'] 
u + by + 



P 



12 



.(59) 



Lq is again the zeroth order Lagrangian giving the action (23) but recast 
in different variables, the tag stands for derivatives with respect to m, Vt is 
some numerical constant, and r is given by: 
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There is also a first order term which vanishes by the EOM, as it must. From 
this we get the EOM for 5y: 



For a plane wave solution with |k| = we finally get: 



Expanding for large u we now find exponential asymptotic behavior with 
one positive exponent, giving the non normalizable mode, and one negative 
giving the normalizable mode. 

Next we solve the equation numerically for the lowest mode exactly as we 
did for a and /3. The result is shown in figure 6. 

In BH embeddings, again there are no normalizable modes. This follows 
from the equation, as expanding it around the horizon (p'^ = 1 + e) we find 
that: 

f 8uj , \ ^ / 8w , \ 
y ^ Acos [ me + 5 sm In e . (63) 



6 Conductivities 

Next we will calculate the conductivities. We will do this in a more general 
case in which there is a B2 field (|9]) . This system is interesting as it might be 
used as a holographic model for the QHE. The reason is that in this model 
we have the following CS term on the D6-brane: 



CiA^sAFAF^ / OFAF. (64) 

J AdSi 

This seems to describe a 3d system of gauge and matter fields with a CS term 



and therefore we expect QHE conductivities analogous to 11 . However this 
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system is supersymmetric and it seems strange that it will have a term such as 



(64). The field theory corresponds to the 3d intersection of the D6-brane with 
various branes preserving A/" = 3 supersymmetry, and so should not generate 
any CS term. This is also apparent from the quiver diagram (Fig. 1). Since 
the added fiavors are a hypermultiplet in the fundamental representation, 
the contributions to the CS term from the fermions should cancel. In order 
to better understand this we will calculate the conductivities. 

Following fT2^, we will add a non-constant gauge field on the probe brane: 

A = { tci + -tijXjb + ai{r) j (65) 

^0 = ^«o(r), (66) 

where Cj, b are the externally applied electric and magnetic fields respectively, 
and Oj, will be given by the equation of motion and they will determine 
the currents. All the fields have been conveniently scaled. 

In order to evaluate the conductivities we need to calculate the action. 
We will start with the DBI part, but before that we need to point out that 
there is also a CS term of the form: 

S = 2n J C-iAB2AF (67) 

that sources the CP^ components of the gauge field for the background of 
interest. So we must also include the terms: 

To simplify matters, I will work with the following ansatz: 

= a^{r) = (a^(r) + C^) cos{d). (69) 

This solves the a,^ EOM, and furthermore it is necessary if one wants ^, oq 
and Oj to be 6 independent. is an arbitrary constant that for simplicity 
I will set to zero. Now I can evaluate the DBI part, which after a lengthly 
calculation gives: 



Sdbi = ^J dr sm2^^1 + {a^ + icos20' ^/fF): (70) 



where 



19 



/(r) = 4 



+ + ((1 

and where we have defined: 



(71) 



and 



/ 



The relevant CS terms are: 



Combining all of these together we get the action: 



S ^ J dr 



sin 2^-^ 1 + {a^ + £cos202a/7(0 



9^6 



(72) 



(73) 



Scs = ^ / ^3 A (^2 + 27rF) A (^2 + 27rF) (74) 
+ (27r)^ / ^ + ^^^^ ^ + 27rF) A (S2 + 2nF). 



+ 



210F 
16A;2 



27^2 --V >./v.v/ 3. 213^2 

{f sin(40e' - av,(a; - 2^sin(20e') - ^cos(20a;) 
cos(2{)(e X a' - 6ao)(a^ + £cos(2^)) 



sin4ee' 
(75) 



As S depends on oq, only through derivatives, there are 3 conserved 
quantities: 



^11 



- ^ 

16k^ dL 
da' 



(76) 
(77) 
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■ dL 



where are the components transverse and parallel to the electric field, 

respectively. It can be shown [12j that these are the physical charge and 
current densities. From this we can, after a prolonged calculation, get the 
following relations: 



^ (p(H--(iW-^))-%fc) ^^^^ 

A//3(r)7(r) - a'^{r) 



^An = jr^ (80) 









e\bp 


p 1 


e 




)-b\e\n 



1 - ■ « 



where: 



«ll W = (81) 



p(r) = p - 6cos2^(a^ + £cos 20 (82) 
(r ) = j_L — I e I cos 2^0^ + i cos 2^ (83) 



7(r) = (i-^)(p2^(^)8i?Sin2 2e(i + 4(a^ + £cos20') 
(j X e)2 + (J.e)2 



e 



2 



».6 



^4) 



a(r)=jxe-6p(l-^). (85) 



^6 
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Inserting this in the DBI part we get: 



R'^ f^^ sin^(20y (r)/3(r) ^^^^ 



where: 



1 - 



H 



1 + 4(a^ + £ cos 20^/1 + (av, + £cos202. (87) 



6.1 BH embeddings 



For BH embeddings the calculation follows as in 12 13 . The argument in 



the square root in the denominator of (87) must not be negative, however 
the positive factor is zero when: 



/3(r)=0^|er = (l-^)(62 + (^f) (88) 

which always has a solution for some > th- The only way the square root 
can be non-negative there is if: 

a(r,) = 0^1, X e-6p,(l- ^) =0 (89) 

^* 

where = p{r^) and = jir^). Furthermore, when r < r^,, (3 becomes 
negative, and the only way the entire square root can stay positive is if 7 
also becomes negative in that range. Therefore, that factor must also vanish 
at that point, so in addition: 

7(r*) = 0^(l-^)-i((j:x)^ + (j.||)^) 

= P*^ + C-^fR^sm^2^,{l + A{a,^ + icos2Q^). (90) 

From this we can derive a relation between the currents and the external 
fields and thus calculate the conductivity. When the dust settles we get: 
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= y^^s (91) 



/ p.^ + )8i?6sin^2e.(l + 4(a^, + f cos2e.)^) b^p,^ 



The conductivity in the BH embeddings is similar to that of a normal metal. 
In particular the transverse conductivity depends linearly on p and ^ (at 
least when b » rjj). 

6.2 MN embedding 

The previous calculation is true only in the BH embeddings, as we must 
assume that the brane fills the all range th < r < oo, otherwise it may end 
before and the argument breaks down. However, in the MN embeddings 
case we have a boundary condition coming from finiteness of the derivatives 
of the components of the gauge field. This is necessary as divergences in 
the derivatives of the gauge field signal that there are sources at the tip. 
This means that there are strings connecting the probe brane with the base 
branes. As it turns out, these strings exert a force on the probe brane forcing 



it to turn to a BH embedding 14 . So in order to have MN embeddings we 
must have no sources, which therefore also requires that the derivatives of 
the components of the gauge field be finite. As we previously mentioned, ^' 



diverges at r = tq. At that point a' must be finite. Equations (79 81) reveal 
that the only way that these terms can be finite is if: 

p(ro) = ^ p = 6 cos 2^{a^ + £ cos 2^) Ij^f = b{i - a^ro)) (93) 
Jll = ^ ail = (94) 

j±{ro) = o-_L = cos2^(a^ + £cos2^)|^^| = i - a^(ro). (95) 

We see that in the MN case we get a vanishing longitudinal conductivity 
and a non-zero transverse conductivity. As MN embeddings possess a mass 
gap we expect them to be insulators, therefore the vanishing of a\\ is ex- 
pected. The non-zero transverse conductivity is surprising and suggests that 
the system is in a quantum hall state. However that is inconsistent with the 
supersymmetry. 
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In order to better understand this we shall look at the EOM, which 



(aV,+^cos2g)sin(20^/(r) 



+ 16 cos(2^) (ba'Q + a' x e) 



First, one can notice that there is the following solution to all the gauge 
field EOM: 



a'o = a' = a^ + i cos 2^ = 0, (97) 

which implies that a± = 0. This solution corresponds to 2ttF = —B2, and 
sets the charge density and currents to zero regardless of e and b. Furthermore 
one can also find MN solutions to the ^ EOM. This solution appears to be 
consistent with what we know about the field theory on the D6-brane. 

However, there may be other solutions with a non-zero value of + 
£cos2^ at the tip, and therefore non-zero currents. These will generate fiux 
on the two-sphere spanned by the coordinates 9, 0, given by: 



27rFe0 + 52,00 = dea^ - — sm{9) cos(2^) 

= -— sin(^)(a^ + £cos(20). (98) 

This is important since this two-sphere shrinks to zero at the tip so the total 
fiux on it must be zero. Again, this is required by the absence of sources, 
which whould render the MN embedding unstable. This then forces + 
£cos2^ = to be zero everywhere, and thus leads to a vanishing transverse 
conductivity. 



7 Conclusion 

We analyzed the ABJM model with fundamental massive matter and showed 
that it exhibits a first order phase transition between solutions which end 
before the horizon (favored at low temperature) and those that enter the 
black hole (favored at high temperature). This is similar to the D3/D7 

•^It is important to note that the equation cannot be derived from the action ( |75[ ), but 
must be derived from the more general Lagrangian, applying the ansatz afterward. 
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system which exhibits the same phenomenon. In the field theory this is a 
transition between an insulator (MN embedding), and a metallic state (BH 
embedding). 

We have shown that the conductivities of the BH embedding are as ex- 
pected in a metallic state. We have also shown that the longitudinal conduc- 
tivity vanishes in the MN embedding as expected. We have further argued 
that the transverse conductivity must also vanish, in accordance with the 
supersymmetry. 
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